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Abstract 

Spectroscopic observations of 51 Pegasi and r Bootis show no periodic changes in the 
shapes of their hne profiles; these results for 51 Peg are in significant conflict with those 
reported by Gray and Hatzes (1997). Our detection limits are small enough to rule out 
nonradial pulsations as the cause of the variability in r Boo, but not in 51 Peg. The absence 
of line shape changes is consistent with these stars' radial velocity variability arising from 
planetary-mass companions. 

Subject headings: line: profiles — stars: individual (51 Peg, r Boo) — stars: oscillations 
— stars: planetary systems 

I. Introduction 

Short-period radial velocity variability of several tens of m s~^ has recently been 
detected in several Sun-like stars, suggesting the presence of planets with roughly jovian 
mass in orbits smaller than 0.25 AU (Mayor & Queloz 1995, Butler et al 1997, Noyes et 
al. 1997). Gray (1997) (henceforth G97) and Gray & Hatzes (1997) (henceforth GH) 
challenged this planetary interpretation, based on evidence for changes, synchronous with 
the radial velocity variation, in the spectrum of 51 Pegasi. The parameters monitored 
by G97 and GH were the shape of the A6253 Fe I line (parameterized by its bisector 
shape), and the ratio of the depths of A6252 V I and A6253 Fe I. Line profile changes 
would not be expected from the gravitational attraction of an orbiting body, but would be 
a natural consequence of dynamic processes (perhaps nonradial pulsations) in the stellar 
atmosphere. It is important to know whether Gray's challenge is correct; if it is, great 
efforts are being expended to explain planets that may not exist. We report here the 
results of an independent spectroscopic study of the stars 51 Peg and r Boo, searching for 
confirming evidence for periodic variations in line shape or depth. In order to get the most 
information out of our moderate-resolution data we have introduced some new analysis 
techniques. These and the related modeling of line profiles that have been distorted by 
nonradial pulsations are described in Brown et al. (1998); here we shall report only our 
results. 

II. Observations and Modeling 

We obtained the observations described here using the Advanced Fiber Optic Echelle 
(AFOE) spectrograph at the 1.5 m Tillinghast telescope at Mt. Hopkins, AZ. This 
spectrograph is designed for precise radial velocity measurements, and employs several 
techniques to assure a stable spectrograph point spread function (Brown et al. 1994). We 
observed 51 Peg and r Boo with the AFOE configured to give \/6\ = R c^ 50000. A typical 
observation of each star consisted of a sequence of 3 sequential integrations, each lasting 
10 m, and each giving S/N of about 150 in the continuum near 600 nm. Because they were 
intended for radial velocity studies, a large majority of our observations were taken using 
an I2 absorption cell. In these data, the stellar spectrum between 500 nm and 610 nm 
is heavily blended with the I2 spectrum, therefore we did not use this wavelength range 
in any of our line shape analysis. In most cases we took only one 3-integration sequence 
each night, but for a few nights we took 2 such sequences, one at the beginning of the 
night and one at the end. The observations of 51 Peg consist of 59 spectra, encompassing 
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20 sequences taken on 18 separate nights between 1995 Nov. 2 and 1996 Jul. 26. Those 
of T Boo consist of 90 spectra, 30 sequences, and 23 nights between 1996 Jun. 25 and 
1997 Mar. 26. The 56 m s~^ radial velocity variation of 51 Peg (period 4.231 d) is clearly 
evident in these data, as is the much larger (468 ni s~^, 3.313 d) variation of r Boo. 

The AFOE's spectral resolution is not high enough to allow useful measurements of 
spectrum line bisectors. We therefore adopted a different method for characterizing the 
shapes of spectrum lines - a decomposition in terms of Hermite functions 7Yi(A/a), defined 
by 

-) = ^exp(-— ^)if„(-; 
a Jo 2(7^ o 



^n(-) = ^exp(-;^)if„(-), (1) 



where Hn{\/(j) is a Hermite polynomial as defined by Abromowitz and Stegun (1972), 
and Nn is a normalization factor given by N^ = (2'^n!y^)~^/^. The decomposition 
involved representing pieces of spectrum /(A) typically 1.5 nm wide and centered on a 
fiducial wavelength A/ as follows: 



/(A) = C[l + ^(A-Aj 
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where C and S describe the brightness and slope of the local continuum, Dj is the depth 
relative to continuum brightness of the j^^ spectrum line, \cj is the center wavelength of 
the j^^ line, and a and the coefficients hi describe the width and the departure from a pure 
Gaussian shape, assumed common to all the lines. The parts of the line profile described 
by hi and /i2 have been absorbed into Ac and a, respectively. Thus, a time series of the hi 
coefficients contains all of the line shape information in our observations, while the ratio of 
the D parameters for the A6252 V I and the A6253 Fe I lines reproduces G97's measurement 
of line depth ratio. In what follows, we express the line shape coefficients hi as a percentage 
of Dj. For the analysis described below, we fit 3 (for 51 Peg) or 4 (for r Boo) sections of 
spectrum, each about 1.5 nm wide. These wavelength ranges contained about 60 distinct 
spectrum lines for each star. This use of multiple lines to gain noise immunity is justified 
when seeking evidence for pulsations, since both theory and observation (of, e.g., 6 Scuti 
stars - see Kennelly et al. 1997) indicate that pulsations modify line shapes in the same 
way for all but very strong spectral lines. 

To compare our results with those of G97 and GH, we were obliged to model the line 
shapes emerging from pulsating and rotating stars, computing both Hermite expansion 
coefficients hi and line bisector shape indices for a wide range of nonradial pulsation 
parameters. The modeling process was similar to that described by, for example, Vogt & 
Penrod (1983), by Hatzes (1996), and by Schrijvers et al. (1997). At each point on a fine 
grid sampling the observable stellar hemisphere, we computed the line-of-sight component 
of velocity resulting from the sum of rotation and pulsation. The line profiles emerging from 
each point were then added together with appropriate Doppler shifts and limb darkening 
to make up the observed profile from the complete star. We then convolved the disk- 
integrated profiles with instrumental point spread functions to give those that would have 
been observed with R = 100000 (to simulate G97's data) and with R = 50000 (to simulate 



the AFOE). Last, we analyzed a series of such hne profiles computed at various phases 
of the oscillation cycle as if they were real data, to determine the relationship between 
measures of bisector shape and our own hi coefficients. 

The principal conclusion from the modeling applies to slowly-rotating stars such as 
51 Peg and to fluid motions (such as those caused by pulsation) that are sensibly constant 
over the range of depths where the spectrum line is formed. In this case, our hs coefficient, 
the bisector span Sb, and bisector curvature Cb all measure essentially the same quantity: 
the lowest-order antisymmetric distortion of the line profile that is not a simple line shift. 
For the residual flux values used by GH in their line bisector studies (.85, .71, .48) and a 
central line depth of 0.6 (appropriate for A6253 Fe I), the proportionality between hs, Sb, 
and Cb turns out to be 

a hs - 0.35Sb ~ -2.2Cb , (3) 

where a is the Gaussian line width parameter in velocity units. For processes {e.g. 
granulation) involving depth-dependent correlations between velocity and temperature, 
this equivalence need not hold. But for the circumstances that most concern us, hs is a 
good proxy for the bisector-based measures of line shape oscillation. We also computed 
the ratio between hs and the RV amplitude. This ratio is quite variable, depending upon 
the mode spherical harmonic degree i, the azimuthal order, the ratio of horizontal to 
vertical velocities, and the inclination of the pulsation axis. It is nevertheless possible to 
characterize the typical behavior for each £, which allowed us to estimate how large an 
hs variation to expect, given a measured radial velocity variation. Finally, we found that 
slowly rotating stars display little change in line shape for a given radial velocity signal, 
but that the shape changes grow more pronounced as v sin i starts to exceed the intrinsic 
line width. Some of these conclusions disagree in important ways with the results of similar 
modeling by GH; for more details, see Brown et al. (1998). 

III. Results of Time-Series Analysis 

A. 51 Pegasi 

We searched for periodic signals in the 51 Peg time series of hi and line depth ratio 
in two ways: we computed periodograms of the time series, and we fit the amplitudes 
and phases of sinusoidal functions with periods of 4.231 d (the radial velocity period) 
and 2.575 d (the period of the largest bisector curvature signal seen by GH). No obvious 
narrow-band signals appear in the periodograms at any frequency, for any of the hi or for 
the line strength ratio. We therefore used the periodograms merely to provide estimates 
of the noise in each of the parameters, using for this purpose the values within a frequency 
range of ±1 cycle per month about the periods of interest. Table 1 shows the results of the 
sinusoid fitting, with error estimates taken from the periodograms as just described. For 
51 Peg, only the h^ amplitudes and the line ratio amplitude for the 2.575 d period exceed 
zero by as much as the estimated error a; the largest difference from zero is 1.7cr. For the 
4.231 d radial velocity period, we measure the variation in line depth ratio to be a little 
less than half that reported by G97, but the uncertainty is such that our result is about 
equally consistent with G97's value or with no variation at all. 

What does the non-detection of an hs signal mean in physical terms? Two kinds of 



answer are possible, as illustrated in Figure 1. The line bisector curvature and the hs 
fitting coefficient measure similar characteristics of spectral line asymmetry. The ratio of 
these two quantities is therefore fairly well defined, and is largely independent of model- 
based assumptions. Thus one can estimate the hs amplitude that corresponds to the 45 
m s~^ curvature amplitude reported by GH. The resulting amplitude is 2.54%, which is 
inconsistent with the measured hs amplitude with lOcr significance. By this measure, our 
observations are therefore strongly inconsistent with those by GH. On the other hand, if one 
assumes a nonradial pulsation with i — 5 and an observed radial velocity amplitude of 56 m 
s~^, then our modeling shows that under typical circumstances the expected amplitudes of 
all hi coefficients should be below 0.5%. Such amplitudes are only marginally detectable in 
our data; therefore, nonradial pulsations cannot be excluded by this test as the cause of the 
radial velocity signal in 51 Peg. The surface fiow velocities associated with this pulsation 
would peak at about 5 km s~^, so the linearity assumption involved in scaling our model 
results would probably be violated. The dissipation associated with such near-sonic fiows 
would likely cause larger line profile variations than those we calculated, however, so the 
actual limit on pulsation amplitude is probably smaller than that given above. 

B. T Bootis 

The star r Boo rotates more rapidly {vsini = 15 km s~^), than does 51 Peg and 
the amplitude of its sinusoidal radial velocity variation (468 m s~^) is a factor of 9 larger 
(Marcy et al. 1997). Both of these circumstances make pulsations easier to detect in r 
Boo than in 51 Peg. If the radial velocity variability of r Boo were to arise from nonradial 
pulsations, then our models show that the anticipated signals in h^, h^, and h^ should all 
have magnitudes of several percent. 

We searched for such periodic signals in the line shapes of r Boo just as we did in 
those of 51 Peg, by examining the periodograms and by fitting sinusoids to the known 
radial velocity period of 3.313 d. The fitted values and their errors may be found in the 
last column of Table 1. Again, none of the fitted sinusoid amplitudes are significantly 
larger than the estimated errors, so no line shape changes at the radial velocity period 
have been detected. Since r Boo is hotter than 51 Peg, its A6252 V I line is too weak to 
allow a meaningful depth measurement. For this reason. Table 1 gives no line depth ratio 
estimate for this star. 

If nonradial pulsations were responsible for r Boo's radial velocity variation, our 
models lead us to expect an hs amplitude of about 15%. Figure 2 shows the observed 
hs values for r Boo phased to the radial velocity period, as well as a comparison between 
the fitted sinusoid and that expected from the pulsation models. In this case the results are 
unambiguous: pulsations with large enough amplitude to cause the radial velocity signal 
would also cause observable distortions of the line shapes, which are not seen. 

The hs periodogram for r Boo shows a wide region of excess power between 0.2 c d~^ 
and 0.4 c d~^. Its amplitude is small, however, (the rms signal associated with this entire 
band is about 0.6%) and the bandwidth is comparable to the estimated rotation frequency 
of the star, which is thought to be identical to the radial velocity frequency. Perhaps this 
hs variation results from magnetic active regions rotating across the disk; it may also be 
associated with the excess radial velocity noise noted in r Boo by Butler et al. (1997). In 



any case, it is too small in amplitude and covers too wide a frequency band to be identified 
directly with the radial velocity signal. 

IV. Discussion 

To summarize our conclusions: (a) We find no evidence for pulsations in 51 Peg. Our 
upper limits on line shape variations are in significant conflict with the bisector curvature 
measurements by GH. (b) For the observed radial velocity signal, our models imply line 
shape variations that are too small to be readily detected, either by our methods or in 
GH's bisector data. Thus, pulsations are not excluded as an explanation for 51 Peg's radial 
velocity variability, (c) The radial velocity variability of r Boo, whose period is similar to 
that of 51 Peg, does not arise from pulsations. 

These conclusions have implications for the study of both extra-solar planets and of 
stellar pulsations. The absence of a pulsation signal in r Boo supports the idea that 
companions of roughly jovian mass may exist in close orbits around Sun-like stars. If one 
such object circles r Boo then it is plausible that another circles 51 Peg, questions about 
the companions' process of origin notwithstanding. As an issue in stellar seismology, 
the absence of large-amplitude (by solar standards) pulsations in Sun-like stars would be 
disappointing but not surprising. Our evidence regarding 51 Peg itself is not conclusive, 
while the observations by GH raise interesting questions but lack conflrmation. Further 
data on this star may therefore be desirable. 

We are grateful to the rest of the AFOE team (Martin Krockenberger and Adam 
Contos) and to the staff at SAO's Whipple Observatory (Bastian van't Saant, Ted Groner, 
Perry Berlind, Jim Peters, and Wayne Peters) for their assistance in obtaining and reducing 
the observations described here. We thank Coen Schrijvers for providing test cases against 
which to compare our line proflle simulation code, and Artie Hatzes for many useful 
discussions. 
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Table 1 

Observed Amplitudes and Errors for 51 Peg & r Boo 





51 Peg 


51 Peg 


r Boo 




4.231 d 


2.575 d 


3.313 d 


h3 (%) 


0.147 ± 0.26 


0.344 ± 0.38 


0.270 ± 0.35 


/l4 (%) 


0.272 ± 0.42 


0.238 ± 0.45 


0.148 ± 0.30 


/15 (%) 


0.219 ± 0.19 


0.323 ± 0.19 


0.190 ± 0.17 


Line Ratio 


0.277 ± 0.30 


0.416 ± 0.28 


- 



Figure Captions 

Figure 1. The line stiape parameter h^ for 51 Peg plotted against phase of the radial 
velocity variation. Filled circles are data points covering one cycle of the RV variation; 
open circles are the same data covering a wider phase range, plotted for clarity of the 
phase relations. The solid curve is a least-squares sine wave fit to the observations, 
with amplitude 0.147%. The other curves are the typical hs signals expected from 
pulsation models yielding (dashed) an RV signal V^op = 56 m s~^, and (dot-dashed) 
a line bisector curvature signal of 45 m s~^. 

Figure 2. Same as Fig. 1, but for r Boo. In this case no bisector curvature information is 
available, so the corresponding (dot-dashed) curve is not shown. 
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